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O pomiarach parametrów funkcji statystycznych stosowanych do opisu szumów 
w obecności tła aparaturowego 


The Determination of the Statistical Parameters of Stochastic Signals 
in the Presence of Apparatus Noise 


Измерение статистических свойств флуктуаций при наличии аппаратурных шумов 


WSTĘP 


Niniejsze opracowanie powstało w wyniku prac nad fluktu- 
acjami /szumami/ prądów termoemisji jonowej w spektrometrach 
паз. 

Strumień jonów powstający w źródle spektrometru ma zwykle 
natężenie mniejsze niż 1078 A. Po wej$ciu и sektorowe pole 
magnetyczne następuje rozdzielenie izotopów i do kolektora 
jonów dociera izotopowy strumień o natężeniu mniejszym niż 
10770 A, Klasyczny sposób rejestracji polega na skierowaniu 
tego strumienia do kolektora jonów i wzmacnianiu we wzmacnia- 


czu elektrometrycznym o dużym oporze wejściowym. Pomiar natę- 
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żenia prądu jonowego odpowiedniego izotopu spośród kilku izoto- 
pów dostarcza wiadomości o składzie izotopowym analizowanej 
próbki, procesach zachodzących na powierzchni enitera, reakcjach 
chemicznych w źródle, wyróżnieniu izotopowym następującym w pro- 
cesie desorpcji termicznej itd. 

Warto w tym miejscu stwierdzić, że z punktu widzenia statys- 
tyki matematycznej, wiadomości te uzyskuje się poprzez pomiar je- 
dynie wartości średniej natężenia ргади /składowej stałej/. W tych 
typowych badaniach pomija się fakt istnienia fluktuacji płynące- 
go pradu jonowego. Zatem słuszne wydaje się mniemanie, że śred- 
nie wyższego rzędu, tzw. momenty centralne rozkładu chwilowych 
wartości anplitud fluktuacji mogą dostarczyć bogatszych wiado- 
mości o zachodzących na powierzchni jonizującej procesach ta- 
kich, jak: desorpcja, dyfuzja powierzchniowa, jonizacja i in- 
nych, zachodzących na powierzchni ciał stałych, w warunkach wy- 
sokich temperatur. Podobnie parametry innych funkcji statystycz-, 
nych /funkcji autokorelacji, funkcji spektralnej gęstości mocy/ 
mogą dostarczyć pełniejszego opisu tych procesów. 


PROBLEMY APARATUROWE 


W przypadku badań fluktuacji prądów jonowych należy sto- 
sować elektrometryczne wzmacniacze szerokopasmowe [1] . Posze- 
rzenie pasma wzmacnianych częstości może być uzyskane przez zas- 
tosowanie wzmacniaczy z ujemnym sprzężenism zwrotnym, W przy- 
padku tzw. 100-procentowego ujemnego sprzężenia zwrotnego, wzma- 
cniacz taki możemy nazwać "przetwornikiem: pręd-napięcie". W ukła- 
dzie takim, dla skzadowej stałej zachodzi: u = WSA zaś dla 


we 
chwilowej wartości fluktuacji: ТУ, = JI), *R,_. Dla wyższych 


częstości wzmacnianego pasma może zachodzić об szef wspói- 
czynnika wzmocnienia wzmacniacza oraz zmniejszenie jego oporu 
wejściowego. Ten fakt należy uwzględnić poprzez odpowiednie ce- 
chowanie toru wzmacniającego i wyznaczenie współczynnika prze- 
twarzania И, /U = И - i/, który maleje monotonicznie wraz ze 
wzrostem częstości. Współczynnik ten, jako funkcję częstości, 
można precyzyjnie wyznaczyć przez modulację strumienia jonów 
sygnałem sinusoidalnym i pomiar powstałej składowej i, SiNWt 
przy różnych częstościach, 
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FUNKCJE STATYSTYCZNE STOSOWANE DO OPISU FLUKTUACJI PRADÓW 


Opis fluktuacji /szumów/ może być dokonany przez podanie 
trzech funkcji statystycznych: spektralnej gęstości mocy, oęs- 
tości prawdopodobieństwa, zwanej inaczej rozkładem chwilowych 
wartości amplitud fluktuacji oraz funkcji autokorelacji. 


Funkcja spekt ra lnej GENS CORE SONS! mo = 
cy S, opisuje zależności częstościowe zawarte v; ѕудпаіе 


szumowym, T 


1 1 
St = lim aka E pom A/ 


Af20 Af Tee T 


Eksperymentalnie funkcję tę wyznacza się, stosując przes- 

trajalny filtr środkowo-przepustowy o małej szerokości pasma 
Af 1 dokonując pomiaru średniej wartości kwadratu fluktu- 

acji dla nastawionej aktualnie częstości бү. Elektroniczne 
uśrednianie następuje po zastosowaniu obwodu RC o dużej sta- 
łej czasowej T. 

Wartość $redniokwadratowa sygnału szumowego w paśmie od 
f, do о jest związana z funkcją spektralnej gęstości mocy w 
następujący sposób: 


f2 

Ge -Í 8, dg- /2/ 
f 
1 


Błąd względny =, popeiniony przy pomiarze Seg zależy od zasto- 
sowanego pasma Óf oraz stałej czasowej T obwodu uśredniające- 
go i 41а przypadku fluktuacji o rozkładzie normalnym wynosi: 


E = МИМАР Т. /3/ 


Podany wyżej błąd względny charakteryzuje część losową błędu. 
№ przypadku widm z ostrą zależnością SĘ od f może wystąpić 
znaczny błąd systematyczny /błąd obciążenia estymatora funkcji 


5/, [1] : 


38 Longin Giadyszewski 


6, = Af? /24 - S"/S, gdzie S" jest drugą 
pochodną funkcji S wzgl;dem zmiennej f. 

w przypadku widm ciągłych, bez ostrych maksimów, część sys- 
tematyczną błędu można zaniedbać. 

W trudnych warunkach eksperymentalnych, gdy poziom mierzo- 
nych szumów jest niewielki, istotne mogą stać się szumy własne 
układu pomiarowego. W takim przypadku można uwzględnić przyczy- 
nek szumu aparaturowego korzystając z twierdzenia : "spektralne 
gęstości mocy szumów nieskorelowanych dodają się”. Zatem, w па- 
szym przypadku: S; = So - Sa’ W ostatnim wzorze Sa oznacza fun- 
kcje S wyznaczoną wprost z pomiaru, natomiast Sa jest spektral- 
ng gestoácia mocy szumu aparaturowego, wyznaczoną w wyniku озо- 
bnego pomiaru, 

Podczas pomiarów fluktuacji termoemisyjnych prądów jonowych 
w spektrometrze mas, okresowo co 20 sek. wyłączano napięcie 
przyspieszajęce jony i układy pomiarowe rejestrowały na taśmie 
samopisu lub za pomocą woltomierza cyfrowego i drukarki poziom 
szumów &pa&raturowych 5„, na przemian z azumami Зе notowanymi 
przy pracujacym źródle jonów. 

w bardzo wielu zjawiskach fizycznych występujące fluktua- 
cje posiadają widmo Srt UN. ЗИ, 3] . Wraz ze wzrostem 
częstości poziom badanych szumów maleje, osięgając przy pewnej 
częstości fe poziom szumów $rutowych. Te ostatnie szumy charak- 
teryzują się szumem białym i dla wszelkich procesów emisji ele- 
ktronów lub jonów dane są wyrażeniem : Ss = 2 e io! fio - nate- 
żenie składowej stałej/. 

Jeśli w badanym zjawisku emisji, szum typu "1/f" i szum 
śrutowy można traktować jako nieskorelowane, wtedy podczas po- 
miarów można uwzględnić również wkład szumu śrutowegoż S56 - Są + 
- Ss . Pomiary takie można prowadzić powyżej granicznej czę- 
stości о aż do momentu, gdy poszukiwana wartość Se osiagnie 
poziom fluktuacji statystycznych przeważającego wtedy szumu 
$rutowego S./ VA? Тт. 

Wielu zjawiskom fizycznym towarzyszą fluktuacje o widmie 
lorentzowskim: 5+ = Зо / 1 + SS. CJ z 2J f. Takie widmo 
obserwuje się w przypadku szumów generacyjno-rekombinacyjnych 
pćłprzewodników, a stała czasowa Т ma charakter średniego cza- 
su życia nośników prądu [4]. 


LE 


O'pomiarach parametrów funkcji... 39 


йа ЕЕ CS autoko re lac j i dostarcza wiadomości 
o zależnościach czasowych, występujących w wielkościach fluktu- 
ujących. Definiuje się ją jako średnię wartość iloczynu "997% 
iy Ту ду • Unormowana funkcja autokorelacji ma postać: 


биде * Суки” ук ati? же: якы /4/ 


Analogowa metoda wyznaczania funkcji autokorelacji polega na 

kolejnym wprowadzaniu wzrastających opóźnień badanego sygnału 

i wyliczaniu średnich wartości iloczynu sygnału nieopóźnionego 

i opóźnionego dla każdej wartości opóźnienia A t..Uzyskany wykres 

funkcji Ry at/ w zależności od At nazywamy autokorelograaen. 
względny błęd standardowy popełniony przy wyznaczaniu ko- 

lejnych wartości funkcji R; at/ można wyliczyć następująco: 


VET \ 2 1/2 
бал М УЗАТ ZUZA ] , 
natomiast dla zerowego opóźnienia: 

Ea T NZ ZPET M /5/ 


Funkcja M t/ procesu stacjonarnego jest funkcją parzys- 
tę: R} at/." LP At/* Jeśli At = O, to dla procesu o wartości 
średniej równej zero funkcja autokorelacji jest równa średnie- 
mu kwadratowi amplitudy szumu бус, = Viae ‚ natomiast dla 
At—» oo wartości chwilowe fluktuacji stają sie niezależne 
1 R, оо" 0.. 

Na podstawie twierdzenia: “funkcja autokorelacji sumy dwu 
sygnałów stacjonarnych i nieskorelowanych jest sumą funkcji 
autokorelacji tych sygnałów” [ 7] , możemy podczas pomiarów 


uwzględnić wkład szumu aparaturowego: R = R, - R,» gdzie 


Zaty 
przez Rz at/ oznaczono funkcję autokorelacji badanego procesu; 
przeż Ra - funkcję autokorelacji uzyskanę na drodze ekspery- 
mentalnej za pomocą aparatury ze znacznymi szumami własnymi 


zaś przez a di funkcję autokoreiacji szumu aparaturowego. 
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Funkcja gęstości prawdopodo- 
bieństwa p może być zdefiniowana jako graniczna war- 
tość stosunku prawdopodobieństwa P wystąpienia fluktuacji o da- 
nej amplitudzie do wielkości przedziału amplitud Ai w którym 
dokonuje się pomiaru. 


onu. РАН 
бөл M T /6/ 


Dla dostatecznie małych wartości “okna analizatora” Ai : 


= P/Aal /7/ 


Pomiaru funkcji gęstości prawdopodobieństwa można dokonać 
przy pomocy analizatorów amplitudy powszechnie stosowanych w 
eksperymentach z dziedziny fizyki jądrowej. Analizatory takie 
dokonują zliczeń dochodzących doń fluktuacji, segregując je we- 
dle skokowo zaprogramowanych przedziałów amplitud /kanałów/. 

Większość sygnałów losowych ma rozkład normalnv. Istnieją 
jednak przypadki sygnałów losowych mających rozkłady znacznie 
odbiegajęce od rozkładu Gaussa |5, 6], dlatego w każdym 
przypadku badania szumów należy przeprowadzić pomiary rozkła- 
dów wartości chwilowych amplitud fluktuacji 1 wykonać test nor- 
malności rozkładu [1] 

Bardzo często uzyskane na drodze eksperymentalnej rozkłady 
amplitudowe nie dają się prosto aproksymować za pomocą znanych 
funkcji rozkładu. Wtedy należy możliwie dokładnie scharaktery- 
zować otrzymany rozkład przez podanie jego tzw. momentów cen- 
tralnych. 

Momentem centralnym k-tego rzędu jp BĘ DZ, EE war= 
tość średnią К-те} potęgi różnicy (н-т)К: к= Ki -n]* » 
gdzie m jest wartością oczekiwaną /średnią/ wyliczang nastepu- 
jeco: m = 1 ЭР}. zaś i, 
i. natomiast p, > prawdopodobieństwo wystąpienia wartości ije 

Dostatecznie dobrą dokładność opisu uzyskuje sie poprzes- 
tając na pierwszych czterech momentach rozkładu eksperymental- 


to możliwe wartości zmiennej losowej 


nego, otrzymanego na drodze skokowego przesuwania okna Л1 ana- 
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lizatora. Wtedy rozkład taki jest histogramem, a zmienna loso- 
wa /szum/ występuje w postaci danych dyskretnych, 

Moment zwykły pierwszy był już definiowany wcześniej jako 
tzw. wartość oczekiwana, Jest on liczony względem początku u- 
kładu bc Roh Md Moment drugi, zwany wariancją Wa 27 

2/4, -n/? p, liczony jest względem wartości średniej m. Moment 
SRM. i czwarty podaje sie zwykle w postaci znormalizowanej 
jako parametry:skośność /ang. eripe s, Sk = / 5 /6? 1 wskaż- 
nik ekscesu /kurtoza/ Ex = A , /6%. 

Dla rozkładów normalnych momenty centralne nieparzyste 
Ар 7O dlan= 1, 3, 5 ..., zaś parzyste: АК = /k-1/G* dla 
k = 2, 4, 6, ..., zatem №. = 62, M , » 36%, itd. Stąd skos- 
ność Sk = O zaś kurtoza Ex » 3, Jeśli dane z eksperymentu stano- 
wią zbiór dyskretnych wartości amplitud szumu, wtedy momenty 
centralne można obliczyć następująco: 


Ак = (1/ м - Vr Gj- mn; 


п) jest tu liczbą amplitud szumu wartości i,, natomiast N jest 


liczbą wszystkich amplitud. a 
Analityczną postać zapisu funkcji rozkładu eksperymentalne- 
go można uzyskać stosując wyrażenia aproksymujące typu funkcji 
Pearsona, uwzględniajęcych cztery momenty centralne EX ign s 
lub wyrażenia w postaci szeregu Edgewortha [13, 14] w któ- 
rym, przy niewielkich odstępstwach rozkładu znalezionego ekspe- 
rynentalnie od rozkładu normalnego : p s 0,5 + 0,5 R; = 1/6 
Sk: ps + 1/24. Ех. P4» gdzie Po jest rozkładen normalnym a Po 
i Pa sa pochodnymi, odpowiednio II 1 III rzędu gęstości normal- 
nej Po’ 
pod przeprowadzania pomiarów rozkładów amplitudowych, 
zmierzających do wyznaczenia momentów centralnych, powstaje bar- 
dzo ważne zagadnienie uwzględnienia wpływu momentów rozkładu 
szumów aparaturowych. Wprowadźmy wielkości Х " zwane semiin `“ 
wariantami [ 10] . wyrażające się poprzez momenty w następują- 
cy sposób; 


Ato Y Af Мј. М 
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ЕЕ 


Jeśli sygneł szumowy będący obiektem badań dodamy do szumu 
aparaturowego, wtedy dla sygnału łącznego 1 z i. * i. możemy 
zapisać wyrażenie, stwierdzające fakt addytywności semiinwarian- 
tów [10] : 


A А. № /9/ 


Powyższe równanie rozpisane dla konkretnych momentów /wzory 8/ 


Ik) 


pozwoli łatwo obliczyć momenty A poszukiwanego rozkładu, gdy 


dysponujemy eksperymentalnie wyznaczonymi momentami Же rozkła- 
du łącznego i momentami | zg rozkładu szumów aparaturowych: 


W. Qm. „0 
e = A с 


/Ao/ 


че (UA. TE 


(a a „27° 
. Mia 0-9 
PU №" Ге 


Jeśli wyznaczymy eksperymentalnie funkcję gęstości prawdopodo- 
bieństwa, to dla każdej snalizowanej amplitudy fluktuacji i. 

w ograniczonym paśmie F = #2 - fi /czyli dla każdego kanału ana- 
lizatora/ możemy wyliczyć błąd względny statystyczny: 


эр 1ДЕ 


W powyższym wzorze T oznacza łączny czas analizy za pomocą ana- 
lizatora wielokanałowego, zaś 41/g względną wielkość prze- 
działu amplitud fluktuacji /wielkość "okna" analizatora/, [111]. 
Błędy popełnione przy wyznaczaniu poszczególnych momentów 
można oszacować według poniższych zależności [8 , 12]: 
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Ga G / VN 
2а „52, 
z -1 . б 
GEN Vztn ] Им 


= V6(N-1 *1)0 N «3 
e V6 IN-1)/(N 1) e37 


2 /11/ 
» a = -1 +3) [N+5 
Ех) = V24N IN 2)N-3) / UN-1)* (N*3) [N+5) 


Twierdzenie Wieners-Chinczyna. 
"Funkcja spektralnej gęstości mocy stacjonarnego procesu loso- 
wego jest transformata Fouriera funkcji autokorelacji tego 
procesu” : 


ж А 
Зы = ^ | zę сонот) от : 
° 
topóźnienie /we wzorze /4/ oznaczona przez At/. Zastosowanie 
transforaaty odwrotnej funkcji S ja daja w razultacie: 


o 
Rer - | 50 cos 0) dt 
0 


W wielu zjawiskach funkcja autokorelacji ma postać wykład 
niczą. Za pomocą relacji Wienera-Chinczyna można łatwo udowod- 
nić, że spektralna gęstość mocy ma wtedy kształt lorentzowski: 


Rq): E2. exp [Ж], $: = 470/11 « (27 uai v је /12/ 


Z twierdzenia Wienera-Chinczyna wynika możliwość wyznacza- 
nia funkcji spektralnej gęstości mocy poprzes funkcję autoko- 
relacji. Operację tę obecnie częściej wykonuje się cyfrowo niż 
analogowo, W pierwszym rzędzie dokonuje się próbkowania anali- 
zowanego przebiegu szumowego. Dane z tej operacji służą do wy- 
znaczenia funkcji autokorelacji. Funkcja ta, w postaci dyskret- 
nych danych, otrzymanych dla różnych, skokowo wprowadzanych o- 
późnień jest wyjściową dla kompucerowego programu transformacji 
Fouriera. 
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Należy zauważyć, że metoda filtracji-bezpośredniego wyzna- 
czania funkcji Sę Jest zupełnie równoważna metodzie funkcji a- 
utokorelacji, pod warunkiem, że proces jest stacjonarny. 

Decyzja о wyborze metody podyktowana jest zwykle przez wa- 
runki aparaturowe i ewentualny dostęp do komputera, 
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SUMMARY 


A short discussion of the analog measurement technique of 
the statistical parameters of stochastic signals in the pre- 
sence of apparatus noise is given. Three main types of sta- 
tistical functions are used to describe the basic properties 
of random signals: probability density function, spectral den- 
sity function and autocorrelation function. 
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РЕЗЮМЕ 


В работе даны замечания об вычислении параметров статис- 
тических функций электрических флюктуаций при наличии аппаратурных 
KyMOB. Описаны: автокорреляционная функция, функция спектраль- 
мой плотности и функция амплитудного расспределения. 

Представлено также дискуссию ошибок совершенных при из- 
мерении этих функций. Представлен способ вычисления вклада аппа- 


ратурных шумов в измеряемые параметры статистических функций. 
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